Higher order intercommutations in Cosmic String Collisions 
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Abstract 

We report the first observation of multiple intercommutation (more than two successive reconnections) of cosmic strings at 
ultra-high collision speeds, and the formation of "kink trains" with up to four closely spaced left- or right-moving kinks. We 
performed a flat space numerical study of abelian Higgs cosmic string intercommutation in the type-11 regime /3 > 1 (where 
/3 = rnl^^i^j./m^auge) up to /3 = 64, the highest value investigated to date. Our results confirm earlier claims that the minimum 
critical speed for double reconnection goes down with increasing P, from ~ 0.98c at /3 = 1 to ~ 0.86c for /3 = 64. Furthermore, 
we observe a qualitative change in the process leading to the second intercommutation: if /3 > 16 it is mediated by a loop 
expanding from the collision point whereas if 1 < /3 < 8 the previously reported "loop" is just an expanding blob of radiation 
which has no topological features and is absorbed by the strings. The multiple reconnections are observed in the loop-mediated, 
deep type-II regime /3 > 16. Triple reconnections appear to be quite generic for collision parameters on the boundary between 
single and double reconnection. For /3 = 16 we observe quadruple events. They result in clustering of small scale structure 
in the form of "kink trains". Our findings suggest that, due to the core interactions, the small scale structure and stochastic 
gravitational wave background of abelian Higgs strings in the strongly type-II regime may be quite different from what would 
be expected from studies of Nambu-Goto strings or of abelian Higgs strings in the « 1 regime. 



Formation of cosmic strings and superstrings is a 
generic outcome of cosmological phase transitions^! and 
some inflationary models, both within Grand Unified 
theories [2J 13] and from branes [U[S]. There is no evidence 
of their existence yet, but strings can lead to a weahh of 
astrophysical phenomena. There is an increasing number 
of surveys and searches looking for observable signatures 
such as CMB anisotropics, gravitational lensing, gravita- 
tional radiation, cosmic rays and gamma ray bursts (see 
classic reviews and more recent updates [5HTU]). 

A string network is expected to reach a scaling solu- 
tion in which statistical properties such as the distance 
between strings become a fixed fraction of the horizon 
size t. The expansion of the Universe pumps energy into 
the string network -by increasing the contribution of long 
strings- but this is balanced by energy losses to radia- 
tion. If these are efficient enough, the contribution from 
the strings to the energy density of the Universe remains 
a small, constant fraction of the dominant form of en- 
ergy (matter or radiation) and is potentially observable. 
Radiation is emitted by oscillating loops, formed when a 
string self-intersects and reconnects, and in bursts. The 
latter are produced by cusps (sections of the string which 
acquire near-luniinal speeds) and, to a lesser extent, by 
kinks created when strings reconnect. The reconnec- 
tion, or intercommutation of strings is therefore an essen- 
tial process that determines and maintains the long-term 
scaling behaviour of the string network. 

Cosmic string intercommutation was first investigated 
numerically in Shellard's pioneering work on global 
strings [H] and later by Matzner[12] for type-II abelian 
Higgs strings with (3 = 2 (using the type-I/type-II termi- 



nology from superconductors, to indicate P < 1 / (3 > 1 
respectively). They pointed out that in ultra-relativistic 
collisions there is a critical center-of-mass velocity Vc 
(collision angle dependent) beyond which strings pass 
through: a loop forms, expanding rapidly from the colli- 
sion point, that catches up with the reconnected strings 
and produces a second intercommutation. A more re- 
cent study [13] focussed on double reconnection, which 
proceeds differently in type-I and type-II strings, with 
the loop only forming in type-II collisions. The angu- 
lar dependence of Vc is calculable from the geometry and 
speeds of the colliding strings [13] . 

In the abelian Higgs model, all evidence to date [T2l[T3] 
indicates that Abrikosov-Nielsen-Olesen (ANO) strings 
with unit winding always reconnect at least once, even 
in ultra-high- speed collisions, with the possible exception 
of near-parallel collisions at very low speed in the type- 
I regime, where the attractive force between the strings 
can make them stick together II 4]. Theoretical predic- 
tions for junction formation based on the Nambu-Goto 
approximation[T5] give an extremely good fit to numeri- 
cal simulations [IB] . 

The simulations in [131 had decreasing resolution with 
increasing /?, and only explored two values in the type-II 
regime (/? = 8,32), but the results suggested that the 
critical velocity for double reconnection would go down 
as a function of /?. This is interesting because one of the 
distinguishing features of cosmic superstrings, as opposed 
to ANO strings, is their low intercommutation probabil- 
ity P ~ 10~^ — 10"-"^ [17], which leads to different scal- 
ing properties, in particular to denser networks (although 
the effect may be weaker than the pt^ ~ 1/P one might 



expect [TS])- However, if the critical velocity of strongly 
type-II ANO strings decreases and becomes comparable 
to the average velocity of the network, the effective inter- 
commutation probability could be much less than one. 

While investigating the critical velocity for large /3 we 
discovered higher order intercommutation events, that is, 
particular collisions in which the strings exchange ends 
more than twice. These we describe now. 

Simulations 

The abelian Higgs model -the relativistic Ginzburg- 
Landau model- involves a complex scalar field (f> and a 
U(l) gauge field with field strength Ff^i, — df_,Au — 
d^Af^ (/i, J/ = 0, 1, 2, 3). In natural units ?i = c = 1 it is 
described by the lagrangian 

The ground state has |(/)| = 77 and zero magnetic field. 
There are two mass scales: scalar excitations have 
mscaiar = V Xt] and gauge field excitations, rugauge = 
\f2er\. Classically, the only relevant parameter in the dy- 
namics is their ratio, ^ = {mscaiar I mgauge)^ = A/2e2, 
which also characterizes the internal structure of the 
ANO vortices EOj . In this paper we are interested 
in the /3 > 1 regime, analogous to a type-II superconduc- 
tor, and in this case the vortices have an inner "scalar" 
core of radius ~ "^Zcaiar which the scalar field departs 
from its vacuum value and vanishes at the center. This is 
surrounded by a larger, "gauge" core of radius ~ 
where the magnetic field is non-zero. As a result, the 
interaction between parallel type-II strings is repulsive. 

We follow the same numerical strategy as in previous 
studies [121 Us] : we use a lattice discretization and place 
a superposition of two oppositely moving ANO strings on 
a three dimensional lattice. This configuration is evolved 
using a leapfrog algorithm. The initial configuration is 
determined by two parameters: the center-of-mass speed 
V of the strings when they are far apart and the angle 
a between them (every collision can be brought to this 
form by an appropriate Lorentz transformation [llj). We 
also impose "freely moving" boundary conditions: after 
each round the fields inside the box are updated using the 
equations of motion, and the fields on the boundaries are 
calculated assuming the strings move unperturbed and 
at constant speeds at the boundaries. All simulations 
shown were done on a 200 x 200 x 400 grid, with lattice 
spacing a = 0.2 and time steps At — 0.02, so the Courant 
condition (here At < holds. 

Our simulations are optimized for the deep type-II 
regime. By solving the two-dimensional, static vortex 
equations one finds, that in a static straight cosmic string 
half of the potential energy in the scalar core is con- 
tained within a radius V^m^^^^^^fdS), where / is a slowly 



varying function with /(I) = l,/(64) = 1.4. Lorentz- 
contraction gives an extra factor j{v)^^ in the direction 
of approach, with j{v) — — v^. This is the small- 

est length scale that has to be resolved. We fix its size 
to match the /3 = 1 case (without loss of generality we 
take A = 2, rj — 1). The scalar core is resolved by at 
least three lattice points up to a center of mass speed of 
V ~ 0.94 — 0.96, which is indicated explicitly in figure [T] 
We define the gauge field core radius to be larger by a 
factor ^/p. The initial string separation is fixed to five 
times the gauge core radius . For /3 < 64 the gauge cores 
still do not overlap in the initial configuration. 

Results 

We simulated the collision of cosmic strings at 
/3 = 1,3.9,4.0,4.1,8,16,31,32,33,49 and 64 for various 
speeds v and angles a to find the critical velocity above 
which the strings effectively pass through each other. 
The results for (3 — 16 are shown in figure [l] other (3 
values and a detailed discussion can be found in ref . pT] . 

Our results confirm a downward trend with /3 of 
the lowest critical velocity Vc.min- Specifically, wc find 
W,vc.rmn) = (4,0.92), (8,0.92), (16,0.90), '(32,0.88), 
(49,0.88) and (64,0.86). 

A closer look at the interaction process reveals that 
the strings do not always intercommute once or twice, 
as previously observed, but also three and four times for 
particular values of the initial speed v and angle a. An 
odd number of reconnections results in overall intercom- 
mutation of the strings segments, and an even number 
in the strings effectively passing through, so we can still 
speak of a critical velocity for the strings passing through. 
However, each reconnection creates small structure on 
the strings in the form of a left- and a right-moving kink. 
In some cases it is not easy to distinguish between one 
and three, or between two and four reconnections, just 
by looking at the energy isosurfaces in the intermediate 
state. But the resulting kinks are clearly visible in the 
final state and can be counted; in case of doubt we use 
this criterion (see fig. [s]). Successive reconnections there- 
fore lead to left- or right-moving "kink trains" , groups of 
up to four closely spaced kinks (one for each intercom- 
mutation). The inter-kink distance within these trains is 
a few string widths, at formation (see figs [2] and |3| . 

The intercommutation process for type II strings un- 
folds as follows. After the collision in which the strings 
exchange ends for the first time two things can happen: 

• For /3 > 16, an expanding loop forms after a short 
delay. U v > Vc, the loop catches up with the two orig- 
inal strings and these reconnect again through the loop. 
This creates a highly curved central region in each string, 
that move towards each other and can sometimes medi- 
ate a third reconnection (see fig. [2]). At this point there 
are two almost anti-aligned string segments and, if they 
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FIG. 1: The number of intercommutations for a range of 
collision velocities and angles for /3 = 16. The symbols A, +,• 
and ★ stand for 1,2,3 or 4 intercommutations respectively. The 
dotted line is a two-parameter fit based on the Nambu-Goto 
approximation[T3|. The dashed line includes core interaction 
effects at high collision angle ~ 180° [21]. Simulations above 
the horizontal line resolve the scalar core size by less than 
three lattice points and are therefore less reliable. Triple and 
quadruple reconnections occur on the boundary betwen single 
and double reconnections. 




FIG. 2; Scalar field isosurfaces with < 0.4. After the 
first intercommutation a loop emerges that catches up with 
the receding strings (second intercommutation), creating two 
highly curved central regions. These move towards each other 
and produce a third intercommutation. Two sets of three 
closely spaced, left- and right-moving kinks are clearly visi- 
ble on each of the strings (indicated by arrows on the upper 
string). (/? = 32, a = 122.7, v = 0.88), t = 3, 5, 10, 17, 20, 22 



are receding sufficiently slowly, a fourth reconnection is 




FIG. 3: Typical string configuration after the fourth inter- 
commutation. The arrows show four closely spaced kinks 
moving up each string; the two corresponding down-moving 
trains are visible below the collision point. (/3 = 16, a = 
122.7, w = 0.92), same isosurfaces, £ = 20. 




FIG. 4: After the first intercommutation a radiation blob 
emerges. The blob catches up with the bridge and is ab- 
sorbed. However, before all radiation is absorbed, a loop 
seems to be formed. This "loop" is not topological, it breaks 
resulting in "Dracula's teeth". (/3 = 4, a = 122.7, i; = 0.98) 
t = 2,4.5,6.5,8. 



possible (see fig. [s]). This is the largest number of inter- 
commutations we have seen, and only for /3 = 16 (while 
triple reconnections are generic). On the other hand, if 
V < Vc the second reconnection does not take place be- 
cause the string loop does not catch up with the strings. 
In this case it will contract again and eventually decay 
into radiation, sometimes after one or a few oscillations. 

• For 1 < /3 < 8 the energy isosurfaces look some- 
what similar, but they reveal a very different intermedi- 
ate state. The "loop" in figure |4] is just a blob of radia- 
tion with no topological features: the (covariant) phase 
of the Higgs field shows no net winding around the appar- 
ent vortex. This is clearly visible in the third timestep, 
where the "loop" meets the string bridges, breaks and is 
absorbed - a real vortex loop would not be able to break. 
In this case the maximum number of reconnections is two. 



Discussion 

The interaction of the cores appears to play an im- 
portant role in interpreting the results. In the type-II 
regime, parallel strings repel. When strings approach 
each other, the interaction between the magnetic cores 
produces a torque that tends to anti-align the string seg- 
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ments, so the actual collision angle is higher than in the 
initial configuration. However, due to this antialignment, 
the collision results in the annihilation of a larger segment 
of string so in fact the strings behave as if they that had 
collided with a smaller angle (which is somewhat coun- 
terintuitive and is explained in [21] )■ For collision angles 
approaching 180° this leads to a different prediction for 
the angular dependence of the critical velocity (fig. [I]) 

Our results appear to confirm the claim [13] that 
the critical speed goes down with increasing /3, from 
Vc ~ 0.96 at /3 = 1 to 0.86 at /3 = 64, ahhough 

this reduction cannot go on indefinitely. A crude (uni- 
versal) lower bound for Vc follows from energy conserva- 
tion: if the strings anti-align locally and a portion L of 
each string is annihilated, the maximum energy available 
to the loop is 2Lj = 2ttR, with R the maximum loop 
radius. If i? < L/2 the second intercommutation cannot 
take place, which happens for v < -y/l — 4/7r2 ~ 0.77. 
These values of Vc are to be contrasted with the average 
root mean square velocity of the network, which has not 
yet been investigated for type-II strings (known values 
range from Vrms — l/\/2 ~ 0.71 for Nambu-Goto strings 
in flat space, to 0.63 — 0.51 in field theory simulations 
with /3 = 1 with cosmic expansion [531 US])- So, dou- 
ble intercommutations may be less rare than in type-I 
collisions, but they are still rare events. However, once 
a kink train is formed its decay time is comparable to 
that of a single kink, and because of its microscopic size 
(it is only a few core widths in length) it is very unlikely 
to be disrupted by intercommutation with another string 
segment. We conclude that the small scale structure of 
strongly type-II abelian Higgs string networks could be 
more clustered than the predictions based on the Nambu- 
Goto approximation with P = 1 [21H^ . The effect on 
cosmological signatures is hard to predict, as multiple 
intercommutations affect energy loss mechanisms in sev- 
eral, competing, ways. First, if a string self-intersects, 
an even number of intercommutations frustrates loop for- 
mation. Second, kinks tend to suppress the formation of 
cusps. Third, kink trains would have a stronger grav- 
itational wave signal than single kinks and would give 
a contribution to the stochastic background not consid- 
ered previously [27H31j . Finally, the contribution from 
reconnections is also enhanced in two ways. While we 
confirm that strongly type-II ANO strings always recon- 
nect {P = 1), an effective intercommutation probability 
P < 1 due to multiple reconnections will still lead to 
denser networks and therefore more reconnections. Most 
of these will be at low velocity, where we expect strong 
anti-alignment, causing longer segments of string to be 
annihilated into radiation than for /? « 1 strings. So we 
would expect stronger and more frequent bursts of radi- 
ation and cosmic rays than for other string types (lower 
/3 and also superstrings) where reconnection bursts are 
always negligible or subdominant |32j . Further work is 
needed to understand these effects quantitatively, which 



affect the cosmological bounds. If cosmic strings are ever 
observed, a crucial question will be whether we can in- 
fer their microscopic structure from the observations, in 
particular whether they are just topological defects or the 
first evidence of superstring theory in the sky. The re- 
sults presented here suggest that, even in the relatively 
simple case of ANO strings, the small scale structure may 
still hold surprises. 
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